Asymmetric quark pairing and color neutrality in the 
Polyakov-Nambu-Jona-Lasinio model of QCD 
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We investigate the effects of confinement and local color neutrality on quark pairing at interme- 
diate chemical potentials in the three flavor Nambu-Jona-Lasinio model with Polyakov loop. While 
prior studies have indicated the presence of a color-flavor-locked (CFL) phase at low temperatures 
and intermediate to high chemical potentials, we find that in the absence of a local color neutrality 
constraint the inclusion of the Polyakov loop gives rise to a new preferred phase in which all quark 
colors and flavors pair, but with unequal magnitudes. We study this asymmetric color-flavor-locked 
(ACFL) phase, which exists even for equal mass quarks, by identifying its location in the phase 
diagram, the order of the associated phase transitions, and its symmetry breaking pattern, which 
proves to be the intersection of the symmetry groups of the 2SC and CFL phases. We also investi- 
gate the effects of the strange quark mass on this new phase and the QCD phase diagram generally. 
Finally, we analyze the effect of a local color neutrality constraint on the QCD phase diagram and 
observe that under at least certain conditions it eliminates both the 2SC and ACFL phases by ren- 
dering them energetically unfavorable. Thus, we observe the previously proposed low temperature 
critical point emerge with the associated continuity between hadronic and color superconducting 
quark matter. 



I. INTRODUCTION 

The phase structure of strongly interacting matter 
has seen an explosion of activity in recent years as the 
boundaries of our experimental probes have continued 
to expand [H-Q. As facilities such as the Large Hadron 
Collider (LHC) and the Relativistic Heavy Ion Collider 
(RHIC) probe matter of ever higher densities and tem- 
peratures, we are able to continually test and refine our 
theoretical models and understanding of matter under 
the extreme conditions encountered in the moments af- 
ter the big bang and in the cores of neutron stars. 

While the fermion sign problem largely restricts the 
techniques of lattice QCD to zero density, one method 
for describing strongly interacting matter throughout the 
phase diagram is the use of effective field theories which 
are built upon the symmetries of QCD. One model which 
has proven useful in this context is the Polyakov-Nambu- 
Jona-Lasinio (PNJL) model, which was developed to de- 
scribe dynamical chiral symmetry breaking and has been 
extended to include quark pairing, confinement, and the 
QCD axial anomaly |6l-[l9|. 

An aspect of the QCD phase diagram of particular 
interest is the nature of quark pairing at intermediate 
chemical potential, /.t. While it is known that for three 
quark flavors a color-flavor-locked (CFL) phase, in which 
all quark flavors and colors pair, is energetically favor- 
able for asymptotically large /i, the preferred pairings for 
/i not asymptotically large are not determined. Calcula- 
tions indicate phases in which only two colors and flavors 
pair (2SC) |2p|. in which one flavor pairs with all others 
(uSC, dSC) [21|, and a phase which has properties of 
both free quarks and hadrons (quarkyonic) [1, IH, . 

In this paper we build on prior studies of the ef- 
fects of confinement on quark pairing in the three fla- 



vor PNJL model by considering a wider range of pairing 
schemes than the CFL and 2SC phases previously con- 
sidered [H, '205. In particular, by permitting distinct ud, 
us, and ds pairing amplitudes, we allow for the possibil- 
ity that the confining mechanism of QCD may not treat 
quark flavors, even for equal masses, on an equal foot- 
ing. Further, by considering a range of strange quark 
masses we investigate the combined effects of this poten- 
tial asymmetry and the decoupling of the strange quark 
sector with increasing mass. 

We also investigate the implications of a local color 
neutrality constraint on the phase structure of dense 
quark matter. While QCD has the capacity to dynami- 
cally achieve local color neutrality by means of a gluon 
field condensate {A'^), the PNJL model lacks the neces- 
sary giuonic degrees of freedom to achieve such neutrality 
in a phase of asymmetric quark pairing (e.g., 2SC, uSC). 
Thus, one must impose such neutrality "by hand" in or- 
der to avoid the large color-electric forces which would 
result from color accumulation (23 - l28| . Prior studies of 
the axial anomaly's influence on the phase structure of 
dense quark matter in the (P)NJL model have either fo- 
cused on pairing structures which are trivially color neu- 
tral P, [13 or have allowed for locally colored phases (2^ . 
By introducing an effectively color-dependent chemical 
potential we impose local color neutrality and study its 
affects on the low temperature portion of the QCD phase 
diagram, most notably its supression of phases of asym- 
metric quark pairing and the subsequent realization of 
quark-hadron continuity. 

The outline of the paper is as follows. We begin in 
Sec. |TT] by recalling the three flavor PNJL model with 
axial anomaly. In Sec. IIIII we construct the phase di- 
agram for massless QCD and identify a new asymmet- 
ric color-flavor-locking (ACFL) phase in which all quarks 
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pair, but with unequal magnitudes. In Sec. IIVI we con- 
struct the phase diagram for massive QCD with various 
strange quark masses in order to study the effects of the 
strange quark mass on both the ACFL phase and the 
phase diagram generally. In Sec.|V]we consider the ACFL 
phase more carefully by studying the associated phase 
transitions and symmetry breaking patterns. Finally, in 
Sec. lVII we impose a color neutrality constraint and inves- 
tigate the resulting suppression of phases of asymmetric 
quark pairing at low temperatures. 



II. THREE-FLAVOR PNJL MODEL 

A. Lagrangian 

The Lagrangian for the three-flavor Nambu-Jona- 
Lasinio model with Polyakov loop at temperature T 
is fliOl 



TABLE I: Coefficients of the Polyakov-loop potential 



C = q{ll 



where the covariant derivative = — if^JJ^o couples 
a static homogeneous gauge field to the quark field 
q and rh is the bare quark mass matrix in flavor space. 

and /C^^-* are effective four- and six-quark interac- 
tions respectively, and (<&,<!>, T) is the Polyakov loop 
potential, which governs the deconfinement transition in 
the pure-gauge sector. 

The four-quark interaction is invariant under the 
SU(3)^ ® SU(3)^ ® U(l)3 ® U(l)^ symmetry group of 
classical QCD, while allowing for spontaneous breaking 
of chiral symmetry and diquark pairing: 

£(4) ^ 8GTr(<?!>V) + 2iJTr(4c^i? + d^dL), (2) 

where ipij = ilRyailLYa is the chiral operator 
and {dR)l = tabc^tjk{qR)iC {qR)^^ and (di)^ = 
£abc£ijk{qL)iC{qL)c are diquark operators of right- and 
left-chirality respectively, with C the charge-conjugation 
operator. The labels a, b, c and «, j, k index color and fla- 
vor respectively. We take G,H > 0, which corresponds 
to attractive four-quark interactions. 

The six-quark interaction reflects the QCD axial 
anomaly by explicitly breaking U(l)^, while retaining in- 
variance under the remaining (physical) QCD symmetry 
group: 

C'^^^ = -SKdet (j) + K'TT[{dldL)(j)] + H.c. (3) 

The final ingredient in our model is the Polyakov loop, 
which serves as an order parameter for confinement in 
the pure-gauge sector [l^ [l^ 



$(x) 



1 



; Tr V exp 




drAoir, x) 



(4) 



where V is the path-ordering operator and (3 = 1/T. 
Writing = A° Aa/2, where the Aq are the Gell-Mann 
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matrices, and working in the Polyakov gauge, in which 
Aq is diagonal, yields = 0Ji^i + 't'&^s- 

As discussed at length in [l|, [ij] , in order to ensure a 
real thermodynamic potential we restrict our attention to 
the case cjfs — 0. Making the standard finite-temperature 
replacements t — > and Aq — > IAq, evaluating Eq. 
(j4|) explicitly for a homogeneous gauge field yields the 
relation 



l + 2cos(/303) 



(5) 



Following Fukushima we describe the pure-gauge decon- 
finement transition via the potential [5, 15] 
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+4($3 + $3)_3($$)2]^ (6) 
where the temperature-dependent coefficients are 

a{T) ^ ao + a, (^^^ + a2 (^^J , b{T) ^ h (^^' ' 

and the and 63 are chosen to correctly reproduce lattice 
QCD resuhs (see Table IJ). In addition, while Tq = 270 
MeV is the critical temperature for the deconfinement 
transition in the pure-gauge sector pTl [l^ , when quarks 
are included in the PNJL model, the transition temper- 
ature deviates from Tq. Therefore, in what follows we 
consider Tq as a parameter of our model, which we will 
set by matching the deconfinement transition atu = 0, 
defined as a maximum in d^/dT (as discussed in pollsol l). 
to the lattice QCD value of T^'^^ = 176 MeV. 



B. Thermodynamic Potential 

Working in mean field, we consider the scalar chiral 
and diquark condensates 



{q^C^zUljq) = diSij. (7) 



Note that there is no sum over z; rather, the right sides of 
Eq. ([7]) are diagonal matrices in flavor space with three 
distinct elements. As shown in P, 0j the mean field 
Lagrangian becomes 

3 

Cmf = ~ ^'^J + ('^ + *'^3A3)7") q] 
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^^A*Aq'c^5tjijq) + ii.c.] (8) 
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where the tj and Ij are the antisymmetric Gell-Mann 
matrices in flavor and color space respectively, V is given 
explicitly below, the effective mass of the jth quark flavor 
is 

K' 

Mj = nij - 4G(Tj + K\ejki\crkcri + — \dj\^, (9) 



and the jth pairing gap is 
A. 



21 H-^a,]d,. 



(10) 



We choose {u,d,s) and {r,g,b) as flavor and color bases 
and use the Gell-Mann matrices X7,5,2 as a representation 
of ii,2,3 and li,2,3, so that Ai represents the dgSb and dbSg 
pairing gap, while A2 and A3 represent the gaps of u^Sft 
and UbSr and Urdg and Ugdr pairs respectively. 

Introducing the Nambu-Gor'kov spinor — 
{q C'(fy/\/2, we may recast our model as a free 
theory with C = '^S~^'^ ~- V — U, where the inverse 
propagator in Nambu-Gor'kov and momentum space is 



s-Hk) = 



-M - ^'7° 



(11) 



with ^' ~ ^ + {4)3X3 and where the sum over j in the off- 
diagonal elements is implied. The condensates directly 
contribute a potential 
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Integrating over the Nambu-Gor'kov fields and per- 
forming the resulting Matsubara sum yields the thermo- 
dynamic potential: 
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ln(l 
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where the En are the 72 poles of the inverse propagator 
in Eq. AE,, = E,, - Etj"" is the difference be- 

tween the eigenvalue and its non-interacting value (with- 
out absolute value), and the factor of 1/2 accounts for 
the double-counting of degrees of freedom in the Nambu- 
Gor'kov formalism. 

Note that in Eq. (IT5|) we have introduced a high- 
momentum cutoff A to regulate the integral. The value of 
A, along with the coupling constants G and K ^ are ini- 
tially fit to empirical mesonic properties and are given 
in Table HI] (parameter set I). Following Abuki et at, 
as we adjust the strange quark mass and coupling K' 
(parameter sets II-IX), rather than recalculating G by 
again fitting the mesonic quantities, we instead, for the 
sake of simplicity, choose G to yield a fixed value for 
{M^ + Ma)/2 = 367.5 MeV The quantitative effects 
of this choice are negligible for the present purposes. 



TABLE II: Parameter sets for the three-flavor PNJL moiiel: 
the strange quark bare mass mg, coupling constants G and K' , 
and Polyakov loop parameter Tq, with a spatial momentum cutoff 
A = 602.3 MeV [9||. Also shown is the constituent strange quark 
mass at /i = T = 0. *In parameter set I all bare quark masses are 
set to zero. In all others we take rriu = 2.5 MeV and = 5.0 
MeV [H. 





ms (MeV) 






To (MeV) 


Ms (MeV) 
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1.926 


12.36 


210 


355.1 


II 


5 


1.928 


12.36 


208 


369.4 


III 


5 


1.928 


51.91 


208 


369.4 


IV 


20 


1.915 


12.36 


207 


392.2 


V 


20 


1.915 


51.91 


207 


392.2 


VI 


40 


1.899 


12.36 


206 


417.5 


VII 


40 


1.899 


51.91 


206 


417.5 


VIII 


80 


1.877 


12.36 


204 


476.6 


IX 


80 


1.877 


51.91 


204 


476.6 



Finally, we consider two values of i^': {!) K' = K , 
which is suggested by applying a Fierz transformation to 
the instanton vertex [17|, and (2) K' = 4:.2K, which al- 
lows the realization of a low T critical point and provides 
for easy comparison to the current literature [ItI . [20j . 
The remaining couplings are HA^ = 1.7 A and KK^ = 
12.36 [ill,[ll[l3. 



III. MASSLESS QCD PHASE DIAGRAM 
A. Without Confinement 

In this section we discuss the phase structure of mass- 
less QCD before moving on to consider the case of three 
different mass quarks. This will allow us to investigate 
both the general effects of quark mass on the phase dia- 
gram, as well its particular influence on a possible ACFL 
phase. Note we do not impose color neutrality in either 
this section or the following, in order that we may con- 
sider the effects of this additional constraint in Sec. IVIl 

When we turn off confinement by setting (j)3 — Q 
and dropping the potential lA{(j)3,T) the thermodynamic 
potential reduces to that considered by Easier and 
Buballa [2^ . The only significant difference between the 
massless NJL model considered here and the massive case 
is that for massless quarks the chiral phase transition is 
first order for all /i, rather than a smooth crossover at 
low ^ (Fig. [iJ . Easier and Euballa have shown that for 
K' > 3.5K a 2SCbec phase appears, which we define as 
a 2SC phase {di — d2 — 0, ds ^ 0) in which Alu^d > fJ-- 
This phase is similarly visible in Fig [TJ separated from 
the hadronic NG phase by a second order phase transi- 
tion, and from the 2SCbcs phase (in which Mu^d < m) by 
a first order transition. 

Anticipating the ACFL phase discussed in Sees. IIVI 
and |Vl in Fig. [5] we show the single diquark conden- 
sate of the CFL phase as a function of temperature 
for fi = 500 MeV. We note that it is roughly constant 
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FIG. 1; (color online) Phase diagram for the NJL model (no 
confinement) with three massless quark flavors. Thick (red) lines 
denote first order transitions while thin (green) lines denote sec- 
ond order transitions. The first order 2SCbec-2SCbcs transition 
is defined by M„ ^(/i, T) = /i. 
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FIG. 3; (color online) Phase diagram for the PNJL model with 
three massless quark flavors. Line types have the same meaning 
as in Fig. [T] with the additional dotted (blue) line representing 
the doconfinement crossover. 
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FIG. 2: (color online) Dimcnsionlcss diquark condensate 
di = di/h!^ in the non-confining NJL model for = 500 MeV. As 
indicated in Fig. [T] the system undergoes a second order phase 
transition from the CFL phase to the QGP at 73 MeV. The lin- 
ear approach to zero for 71 MeV < T < 73 MeV is due to an 
effective cr|dp coupling [ll. |32|. [3^. 



for T < 30 MeV, and then falls as d - ^/T~T for 
30 MeV < T < 71 MeV, before finally vanishing as 
d ^ {Tc ~ T), due to the effective cr|dp coupling induced 
by the axial anomaly [H, HI, . 



parameter sets used in Table |TT1 Minimizing f2 with re- 
spect to the condensates and Polyakov loop, we obtain 
the phase diagram shown in Fig. |3l As has been widely 
reported, the inclusion of the Polyakov loop pulls the chi- 
ral transition to higher temperatures (from 151 MeV to 
193 MeV) , si gnif icantly enlarging the region of symmetry 
breaking [ul, 111, Hi • 

One important consequence of the increase of T^'^^ is 
that the Polyakov loop gives rise to a much larger region 
of 2SCbcs, which we define as a 2SC phase [di = ^2 = 0, 
da ^ 0) in which M^.d < In particular, this phase now 
persists to much higher /i than in the NJL model, where 
it is constrained to roughly 270 MeV < ^ 350 MeV. 

Figure |4] shows the two distinct diquark condensates 
di = d2 and dg for /i = 500 MeV. We find that for T < 20 
MeV, the results arc not significantly altered from the 
NJL model. However, for T > 20 MeV we find that di = 
d2 falls with increasing T, while d^ increases until the 
system undergoes a second order phase transition to the 
2SCbcs phase at 70 MeV, shghtly below the location of 
the CFL-QGP transition in the absence of confinement. 
Thus, the ground state of the system at intermediate 
/i is no longer a symmetric CFL phase, but rather an 
asymmetric CFL phase characterized by < di = d2 < 
da'. 



IV. REALISTIC MASS QCD PHASE DIAGRAM 



B. With Confinement 

In order to construct the phase diagram in the pres- 
ence of the Polyakov loop, we first fix Tq by matching 
the model's deconfinement transition at /i = to the lat- 
tice value of T^^P = 176 MeV. The resulting value of 
To varies slightly with ms, and is given for the various 



Having observed the emergence of an ACFL phase in 
massless QCD, we now consider the effects of realistic 
bare quark masses on both this phase and the phase di- 
agram generally. To do so we construct phase diagrams 
for nis = 0, 20, 40, and 80 MeV. In all cases we take 
rriu = 2.5 MeV and rud — 5.0 MeV, while the coupling 
G is adjusted in order to maintain (M„ -|- Md)/^ = 367.5 
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FIG. 4: (color online) Dimensionless diquark condensates in the 
massless PNJL model for ^ = 500 MeV. 



mass quarks much of the ACFL-2SC phase boundary 
is actually two distinct, but very closely spaced phase 
boundaries. The first boundary, at slightly lower tem- 
perature, separates the ACFL phase from a sliver of a 
uSC phase in which up/down and up/strange quarks 
pair, but down/strange quarks do not. Thus, cross- 
ing this phase boundary corresponds to breaking the 
down/strange quark Cooper pairs. The second bound- 
ary separates the uSC phase from the 2SC and corre- 
sponds to the breaking of the up/strange quark pairs. 
Figure [7] shows these two distinct transitions for exag- 
gerated up and down quark masses (m„ = 0, rud = 40 
MeV, TOs — 80 MeV), in order to make the distinct phase 
boundaries visible. 



V. ASYMMETRIC CFL (ACFL) PHASE 
A. Quark Pairing Amplitudes 



MeV at M = r 0. 

As shown in Figs. [S]and[Sl as the strange quark mass 
increases, the region of ACFL moves to higher /z, effec- 
tively decoupling the strange quark from the up/down 
sector. This is due to the fact that in the limit nia — >■ oo, 
there is insufficient energy to generate strange quarks and 
we are left with an effectively two flavor system. We also 
note that while for small m^, the deconfinement transi- 
tion at large /j, essentially coincides with the breaking of 
up/down quark Cooper pairs (the 2SC-QGP transition), 
as rUs increases the deconfinement temperature moves 
down somewhat. 

Also noteworthy is the fact that except for small TOs 
and K' , a critical point appears on the ACFL-2SC phase 
boundary, separating a first-order transition at lower /i 
from a second-order transition at higher /z. One can sum- 
marize the situation by noting that the ACFL-2SC tran- 
sition is first order when Tc < 50 MeV, and second order 
when Tc > 50 MeV. Thus, for example, for ms = 5, 20 
MeV and K' = K, the phase boundary never drops be- 
low T « 50 MeV and the transition is always second 
order. We note, however, that while the phase boundary 
has a negative slope for large fi, the transition does not 
again become first order when the boundary drops below 
T w 50 MeV. 

As shown by Abuki et al. in the non-confining NJL 
model and by the authors for the massless PNJL model, 
we find that for K' > A.2K, a low T critical point 
emerges [H, [T^- Also, as shown by Easier and Buballa, 
when one allows for 2SC pairing this critical point acts 
as the termination of a line of first order BEC-BCS tran- 
sitions, above which a smooth crossover develops (20l |. 
Interestingly, as shown in Fig. |6l when the 2SCbec phase 
exists, we find that for = 5, 20, and 40 MeV the BEC- 
BCS transition is first order at zero temperature, while 
for nis ~ 80 MeV the critical point drops below the T- 
axis and one obtains a smooth BEC-BCS crossover. 

Finally, while not visible in Figs. [S] and IHl for unequal 



The evolution of color superconducting quark matter 
with increasing temperature can be inferred from Fig. |8l 
At low temperatures, the ACFL phase is essentially iden- 
tical with the CFL phase, with di = 6,2 ~ d^, and has a 
thermodynamic potential well below the QGP. At high 
temperatures, the ACFL phase morphs continuously into 
the 2SCbcs phase, with di = d2 = 0, via a second order 
phase transition. In between these two limiting cases, 
for 20 MeV < T < 70 MeV, the ACFL phase is distinct 
from both the 2SC and QGP phases, and has a thermo- 
dynamic potential below both. 

We also note that while not clearly visible in Fig. |4j 
our calculations indicate that for T > 6 MeV it is always 
energetically favorable to adopt unequal pairing ampli- 
tudes. Thus, while we cannot exclude the possibility of 
a low temperature CFL- ACFL phase transition, it seems 
very likely that the unequal pairing amplitudes persist to 
arbitrarily low temperatures, and that a symmetric CFL 
phase at intermediate ^ is restricted to T = 0. 

We can understand the asymmetric behavior of the 
quark pairing by noting that in our chosen gauge (and 
with 08 = 0) the quark-Polyakov loop coupling is of the 
form 

qA^fq = 4)z{rfr ~ gfg), (14) 

where we have written the color indices explicitly, so that 
the Polyakov loop couples only to red and green quarks. 
Thus, the condensates di (which involves green and blue 
quarks) and d2 (red and blue) are only singly-coupled to 
the Polyakov loop, while (green and red) is doubly- 
coupled. 

One may inquire whether this phase of unequal quark 
pairing is simply an artifact of our choice of — 0, or 
whether such a phase might actually be realized in QCD. 
Unfortunately, in the present model, allowing ^ ren- 
ders the thermodynamic potential complex so that its 
minimization is no longer a well-posed problem. Never- 
theless, our results do demonstrate the possibility of ob- 
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FIG. 5: (color online) Phase diagrams for the PNJL model with realistic up and down quark masses and various 
strange quark masses, where the axial anomaly couplings are taken to be equal {K' = K). Line types have the 
same meaning as in Fig. |31 



taining a phase characterized by non-equal quark pairing, 
and present a challenge to other models of dense quark 
matter to address the question of its realization. 

B. Symmetry Breaking Pattern 

Having identified the region of the phase diagram oc- 
cupied by the ACFL phase as well as the order of the 
associated phase transitions, we next study the symme- 
try breaking pattern of this phase. We begin by noting 
that the symmetry groups of the 2SC and CFL states 
are [ll,[3f 

2SC : SU(2)^g ® SU(2)^ ® SU(2)^ ® U(l)g ® UCl)^, 
CFL:SU(3)^+^+^®Z2, 

where SU(2)rg denotes a rotation in the color subspace 
of red and green quarks, U(l)^ is a rotated baryon con- 
serving symmetry with conserved quantity 

B^Q + h with Q^Q - J—Xs, (15) 



where Is is the isospin operator, Q and Q are the stan- 
dard and rotated (conserved) electromagnetic charge op- 
erators in the 2SC phase, and U(l)s corresponds to mul- 
tiplying the strange quark by an arbitrary phase. 

Since both the CFL and 2SC phases are special cases 
of the ACFL phase, the symmetry group of the ACFL 
phase must be a subset of the symmetry groups of these 
respective phases. Thus, the color-flavor-locking aspect 
of the CFL phase requires that there be no unbroken 
independent color or chiral rotations in the ACFL phase, 
while the SU(2)rg symmetry of the 2SC phase requires 
that there be no unbroken symmetry which mixes blue 
quarks with red or green quarks. A direct calculation 
demonstrates that none of the remaining symmetries are 
broken and we are left with the symmetry group 

ACFL:SU(2)^^^^^^^®Z2. 

In fact, the symmetry group of the ACFL phase is 
simply the intersection of the symmetry groups of the 
2SC and CFL phases. Moreover, this symmetry group is 
identical to that of the CFL phase with unequal strange 
quark mass (s^. Finally, we note that we expect 14 Gold- 
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FIG. 6; (color online) Phase diagrams for the PNJL model with realistic up and down quark masses and various 
strange quark masses, where K' = A.2K. Line types have the same meaning as in Fig. |31 



stone bosons in the ACFL phase, which follows from the 
Sfl + 8l + 1_b — 17 generators of the Lagrangian and the 
3 generators of the ACFL symmetry group. 



VI. COLOR NEUTRALITY 

In the prior sections we have constructed the phase di- 
agram of the PNJL model for both massless and massive 
quarks and have observed the emergence of a new ACFL 
phase at large fji. If our model is to accurately reflect the 
behavior of dense QCD, however, for the homogeneous 
phases which we consider here we must also investigate 
the effects of the requirement of local color neutrality. In 
fact, both the 2SC phase previously reported by Easier 
and Buballa [l^l and the new ACFL phase posess non- 
zero color densities which would, if left unchecked, induce 
large color-electric forces in the superconducting quark 
matter. 

The origin of the net color density, in both the 2SC and 
ACFL phases, is the modification of the quark dispersion 
relations which results from unequal pairing amplitudes 
for red and green quarks compared with blue quarks. In 



the 2SC phase, for example, at fixed particle number the 
pairing of red and green quarks results in a decrease in 
the Fermi energy of these colors. In a system at fixed 
quark chemical potential /i, this results in an increase 
in the density of red and green quarks compared to the 
unpaired blue quarks, and a corresponding net anti-blue 
color density. In QCD, this quark color density is ex- 
actly cancelled by the development of a non-zero expec- 
tation value of the gluon field (i.e., tadpole diagrams), 
and so the homogeneous 2SC phase remains color neu- 
tral [H, |2^. However, having replaced the local SU(3) 
color symmetry of QCD with the global symmetry of the 
PNJL model we now lack the means for dynamically re- 
alizing a neutral ground state. 

The standard method for imposing color neutrality in 
the NJL model is to introduce a set of color chemical 
potential s /ig, w hich are chosen to ensure vanishing color 
densities |26l - [28| : 

= (gtr.g) ^-|^ =0, (16) 

where = Aa/2. In light of our prior discussion, we see 
that the equilibrium value of /ig (i.e., the value required 
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as well as the stability condition 
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FIG. 7; (color online) Partial phase diagram of the three flavor 
PNJL model with K' = 4.2K and (bare) masses of rriu = 0, 
niii = 40 MeV, and rris = 80 MeV. Line types have the same 
meaning as in Fig. [T] 
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FIG. 8: (color onUne) Minimum H/A* vs. T at = 500 MeV 
for the phases: CFL (dot-dash; green), ACFL (dotted; blue), and 
QGP (solid; red) in the massless PNJL model. 



to achieve color neutrality) is proportional to {AD in 
QCD. In both the 2SC and ACFL phases red and green 
quarks pair symmetrically, so we need only include /is, in 
order to ensure that ns = rir + ng — 2nb = 2{nr — n-t) = 0. 
Thus, we modify the Lagrangian from Eq. ([Ij to 



£ = q{il/) — m + /i7" 
-W($,¥,T). 



(6) 



(17) 



In order to obtain the locally color neutral phase di- 
agram we now minimize the thermodynamic potential 
with respect to the condensates ai and d^, and Polyakov 
loop variable (/>3 as before, while imposing the additional 
neutrality constraint 



0, 



(18) 



dns 
dfis 



dfil 



> 0. 



(19) 



Thus, our solution is a saddle point of f2, minimized with 
respect to the condensates, and maximized with respect 
to /is. 

Due to the computational intensity of the saddle point 
problem for our eight variable thermodynamic potential 
we defer a complete assessment of the effects of color neu- 
trality, together with the strange quark mass and con- 
finement, to a future publication. However, we report 
three important results from the massless quark limit at 
T — 0, which give insight into the structure of the full 
color neutral QCD phase diagram. 

First, in the massless quark limit the color neutrality 
constraint eliminates the 2SC phase from a large portion 
of the phase diagram, in favor of a CFL phase. We can 
understand this effect by considering Fig. [Q] At /i = 275 
MeV the thermodynamic potentials of the color neutral 
NG and CFL phases (/is = 0) are nearly equal, indicating 
the location of a phase transition between the NG phase, 
which exists at low fi, and the CFL phase which exists at 
high fi. As the system moves to higher density the energy 
of the NG phase is essentially constant, while both the 
2SC and CFL phases decrease in energy, becoming more 
favorable. 

The crucial affect of the color neutrality constraint is 
visible in the thermodynamic potentials at /t = 285 MeV. 
In the absence of a color neutrality constraint (/*§ = 0), 
we find that the 2SC phase is indeed the lowest energy, 
and therefore the preferred, phase of the system. How- 
ever, in imposing color neutrality, we require the 2SC 
phase to take on a nonzero ns ~ —40 MeV, which results 
in a (physical) 2SC state which is formally higher in en- 
ergy than the colored state. This "additional" energy 
is sufficient to raise fl2SC above both Ung and flcFL, 
with the lower energy CFL phase being the color neutral 
ground state. As the system moves to yet higher fj,, both 
the 2SC and CFL phases continue to move to lower ener- 
gies, with the latter always maintaining a slight energetic 
advantage. Thus, at T = color neutrality eliminates the 
2SC phase altogether. 

A second effect of the local color neutrality constraint 
is the elimination of the ACFL phase at high /t, in favor of 
a symmetric CFL phase. This is somewhat encouraging, 
given our expectation of a CFL phase at asymptotically 
high fi, due to general considerations [13, HI]. Thus, we 
find that the color neutrality constraint disfavors both 
the asymmetric 2SC and ACFL phases, and at least for 
some parameters, leads to the complete exclusion of these 
phases. 

A third important effect of color neutrality, which is 
a corollary of the suppression of the 2SC and ACFL 
phases is the "re-emergence" of a low temperature critical 
point Ail. As shown in |2(|, when one allows for 2SC 
quark pairing (rather than simply a CFL structure) in the 
absence of a local color neutrality constraint this critical 
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FIG. 9; (color online) Minimum thermodynamic potentials in the NG, 2SC, and CFL phases as a function of /ig at 
T = for /i = 275 and 285 MeV, in the three flavor PNJL model with massless quarks and K' = 4:.2K. 



point is eliminated in favor of a second order NG-2SC 
phase transition at intermediate fj,. However, with the 
2SC and ACFL phases eliminated by the color neutrality 
constraint, the system once again realizes quark-hadron 
continuity via a smooth crossover between the NG and 
CFL phases at low temperatures. 

A number of outstanding questions exist regarding the 
PNJL model and the QCD phase diagram which we will 
address in a future publication [39*1 . Foremost among 
them is a complete construction of the QCD phase dia- 
gram which incorporates local color neutrality along with 
realistic quark masses. To the extent that color neutral- 
ity eliminates the 2SC phase from the phase diagram, 
prior results which did not account for 2SC pairing such 
as those in [l], [13] will remain valid. However, the 2SC 
phase will likely persist in smaller regions of the phase 



diagram, at least for some parameter sets; the location 
of these regions and their associated phase transitions 
must be considered carefully. Also the effects of charge 
neutrality and /3-equilibrium, which are important in the 
study of stable quark matter at low temperatures in neu- 
tron stars, have not been completely elucidated in the 
context of the PNJL model. 



VII. ACKNOWLEDGEMENTS 

This research was supported in part by NSF Grant 
PHY09-69790. The authors thank Professor Tetsuo Hat- 
suda for his insightful comments and discussion. Author 
PP thanks Tomoki Ozawa for helpful discussion regard- 
ing symmetry breaking patterns. 



[1] P.D. Powell and G. Baym, Phys. Rev. D 85 074003 
(2012). 

[2] C.G. Callan, R.F. Dashen, and D.J. Gross, Phys. Lett. 

B 78, 307 (1978). 
[3] J. Kuti, B Lukacs, J. Polonyi, and K. Szlachanyi, Phys. 

Lett. B 75 (1980). 
[4] R.D. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 

(1984). 

[5] K. Fukushima and T. Hatsuda, Rep. Prog. Phys. 74, 
014001 (2011). 

[6] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 
(1961). 

[7] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 124, 246 
(1961). 

[8] T. Hatsuda and T. Kunihiro, Phys. Rep. 247, 221 (1994). 

[9] M. Buballa, Phys. Rep. 407, 205 (2005). 
[10] A. M. Polyakov, Phys. Lett. 72B 477 (1978). 
[11] C. Ratti, M.A. Thaler, and W. Weise, Phys. Rev. D 73, 
014019 (2006). 



[12] C. Sasaki, B. Friman, and K. Redlich, Phys. Rev. D 75, 
074013 (2007). 

[13] S. Rofiner, C. Ratti, and W. Weise, Phys. Rev. D 75, 

034007 (2007). 
[14] H. Abuki, Prog. Theo. Phys. Supp. 174, 66 (2008). 
[15] K. Fukushima, Phys. Lett. B 591, 277 (2004). 
[16] M. Alford, Ann. Rev. Nucl. Part. Sci. 51, 131 (2001). 
[17] H. Abuki, G. Baym, T. Hatsuda, and N. Yamamoto, 

Phys. Rev. D 81, 125010 (2010). 
[18] G. 't Hooft, Phys. Rev. Lett. 37, 8 (1976); Phys. Rev. D 

14, 3432 (1976). 
[19] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 44, 

1422 (1970). 

[20] H. Basler and M. Buballa, Phys. Rev. D 82, 094004 
(2010). 

[21] H. Abuki and T. Kunihiro, Nucl. Phys. A 768, 118 

(2006) . 

[22] L. McLerran and R. D. Pisarski, Nucl. Phys. A 796, 83 

(2007) . 



10 



[23] Y. Hidaka, L. D. McLerran, and R. D. Pisarski, Nucl. 

Phys. A 808, 117 (2008). 
[24] A. Gerhold and A. Rebhan, Phys. Rev. D 68, 011502 

(2003). 

[25] D. D. Dietrich and D. H. Rischke, Prog. Part. Nucl. Phys. 

53, 305 (2004). 
[26] K. lida and G. Baym, Phys. Rev. D 63, 074018 (2001). 
[27] A. W. Steiner, S. Reddy, and M. Prakash, Phys. Rev. D 

66 094007 (2002). 
[28] M. Buballa and I. A. Shovkovy, Phys. Rev. D 72, 097501 

(2005). 

[29] Y. Aoki et al, Phys. Lett. B 643, 46 (2006). 

[30] Y. Aoki et al, J. High Energy Phys. 6, 88 (2009). 

[31] J. Beringer et al. (Particle Data Group), Phys. Rev. D 

86, 010001 (2012). 
[32] T. Hatsuda, M. Tachibana, N. Yamamoto, and G. Baym, 

Phys. Rev. Lett. 97, 122001 (2006). 
[33] N. Yamamoto, M. Tachibana, T. Hatsuda, and G. Baym, 

Phys. Rev. D 76, 074001 (2007). 
[34] Note that in our prior work [l[, we did not observe such 

an increase in the region of chiral symmetry breaking. 



This was due to our defining the deconfinement transi- 
tion as the (sUght) discontinuity in d^/dT, which is co- 
incident with the chiral phase transition, rather than its 
maximum finite value. Here, we use the latter definition 
as it more naturally extends to the case of non-zero quark 
masses in which d^/dT is continuous and the deconfine- 
ment "transition" is unquestionably a smooth crossover. 

[35] M. Alford, K. Rajagopal, T. Schaefer, and A. Schmitt, 
Rev. Mod. Phys. 80, 1455 (2008). 

[36] While one often refers to a CFL phase in the case of 
unequal strange quark mass, it is implicitly assumed that 
fj, S> rUs so that di = d2 ~ ds. If M is not much larger 
than TUs then the quark pairing amplitudes will no longer 
be approximately equal and the CFL phase gives way to 
the ACFL phase. 

[37] M. Srednicki and L. Susskind, Nucl. Phys. B 187 93, 
(1981). 

[38] M. Alford, K. Rajagopal, and F. Wilczek, Nucl. Phys. B 

537 443, (1999). 
[39] P.D. Powell and G. Baym (to be published). 



